APS/123-QED 



■ 

o ■ 
o : 

a: 

m ■ 

o\ '■ 
(N ; 

C/2 ! 

73 : 

O 1 

C0 ' 
>V 

^ : 

Ph. 



> 

OO 

o 
o 



o 



X 



THE METHOD OF THE DESCRIPTION OF DYNAMICS NONEQUILIBRIUM 
SYSTEMS WITHIN THE FRAMES OF THE CLASSICAL MECHANICS 

V.M. Somsiko\0 

Laboratory of Physics of the geoheliocosmic relation, Institute of Ionosphere, Almaty, Kazakstan. 

(Dated: February 2, 2008) 

Within the frames of the analytical mechanics the method of the description of dynamics of 
nonequilibrium systems of potentially interacting elements is develops. The method is based on an 
opportunity of representation of nonequilibrium system by set of interacting equilibrium subsystems. 
The equation of motion of interacting subsystems is found. Based on it the Lagrange, Hamilton and 
Liouville equations for subsystems are obtained. The expression of an entropy production is found. 
The way of a substantiation of thermodynamics in the frames of classical mechanic is offered. 
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I. INTRODUCTION 

The investigations of the open nonequilibrium systems 
collided with the big difficulties. One of them is contra- 
diction between classical mechanics and thermodynam- 
ics. Since Boltzmann and up to now the attempts to 
remove this contradiction are undertaken [1-4]. But, as 
a rule, offered solution cannot overcome Poincare's the- 
orem of the recurrence which forbid an establishment of 
equilibrium in Hamilton's systems [2, 3]. The contempo- 
rary explanation of the mechanism of an establishment 
of equilibrium in Hamilton's systems basing on the con- 
dition of mixing is not rigorous enough as it demands 
postulation "coarse grain", i.e. averaging of the phase 
space on physically small volume. The nature of such 
average is impossible to explain in the frame of the clas- 
sical mechanics. 

With the purpose of finding-out of the physical rea- 
son of the contradiction between thermodynamics and 
classical mechanics, we had been analyzed a hard-disks 
system [5, 6]. The non-equilibrium system of disks was 
represented by a set of interacting equilibrium subsys- 
tems (IES). Their analysis has shown that the process 
of an establishment of equilibrium is caused by transfor- 
mation of kinetic energy of IES motion to their internal 
energy in a result of a chaos increasing of the disks veloc- 
ities. I.e. the dissipation for a disks system is determined 
by the transformation of energy of relative motion of IES 
into internal energy by the work of the collective force 
of interaction IES. Therefore the problem about evolu- 
tion of nonequilibrium systems is reduced to determin- 
ing of the forces between IES. For a hard-disks system 
of these forces are easily enough determined due to the 
energy of disks has only a kinetic component. But fun- 
damental forces of interactions is potential, therefore the 
establishment of equilibrium is determined by the forces 
which work redistributes between IES both kinetic, and 
potential energy. Therefore to determine the process of 
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evolution in system, it is necessary to determine, what 
character energy streams between IES created by the in- 
trasystem forces. It is impossible to do basing on the 
Hamilton formalism because formalism is applicable only 
for system as a whole. However streams of energy be- 
tween IES can be determined with the help of equation 
of the systems (EQS) which determine energy exchange 
between IES [7]. 

The EQS allows offering the new approach to the anal- 
ysis of nonequilibrium systems. This approach is based 
on the following assumptions and conditions: 1). The 
closed nonequilibrium system can be submitted by a set 
of motioning in relative to each other of interacting IES; 
2). Dynamics of nonequilibrium system is equivalent of 
the dynamics of interacting IES; 3). The energy of IES 
should be submitted by the sum of internal energy and 
energy of IES motion as the whole; 4). Each element of 
system should be fixed for corresponding IES without de- 
pendence from its mixing in space; 5). During all process 
of a subsystem are considered as equilibrium. 

The applicability of the first condition is guaranteed by 
that it was successfully used earlier, for example, in case 
of a substantiation of a principle of entropy maximum for 
the equilibrium systems [8] . The second condition follows 
from the first. The third condition is necessary for intro- 
duction in the description of dynamics of systems of inter- 
nal energy, as one of key parameters describing dynamics 
of interacting IES. The fourth condition allows avoiding 
a problem of redefinition IES due to particles mixing. 
Last condition is known from statistical physics. It, not 
depriving a task of a generality, removes the problems 
connected to complication of the analysis of redistribu- 
tion of various types of energy between IES. The offered 
approach allows to connect thermodynamics with classi- 
cal mechanics and to analyze some laws of evolution of 
nonequilibrium system. 

Here the substantiation of this approach is offered. The 
EQS is obtained. The expression for the dissipative force 
determining change of internal energy IES is obtained 
also. The analysis of Lagrange, Hamilton and Liouville 
equations for IES is given. Formulas for the entropy and 
it production are submitted. The connection EQS and 
the basic equation of thermodynamics is shown. 
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II. THE EQUATION OF SYSTEM DYNAMICS 

The analytical mechanics as usually constructed from 
a postulate: "virtual work of forces of reaction is al- 
ways equal to zero on any virtual displacement which is 
not breaking set kinematics connections" [9]. Basing on 
this postulate and a condition of monogenic of the active 
forces, come to a principle of the least action, functions 
Lagrange and Hamilton. It is possible to build mechan- 
ics based on the principle of the least action [10]. But 
it is possible to do basing on the law of conservation of 
energy without the requirement monogenic forces. Here 
by this way the equation of the systems interaction will 
be obtained. It will help us discover that the dynamics of 
the system is determined by two types of energy: energy 
of motion of system as the whole and internal energy. To 
each of this type of energy the force which change this 
type are corresponds. 

Let us obtain the dynamics equations for elements and 
their systems in a field of external forces. If time is ho- 
mogeneous, for energy of system, E, we have: E = 0. 
In non-homogeneous space at presence of interaction of 
particles, energy of system depends on both velocities of 
elements, v, and from their coordinates, r. In this case 
the equality, E = 0, takes place only when the energy de- 
pends on two additive parts. One part should be a func- 
tion of velocity, and the second one - coordinates. I.e. the 
energy can be written down as: E — <p{T + U] = const, 

N 

where T = ^Ti(v 2 ), vt is a velocity of i- element, T is 

i=l 

the kinetic energy of the system, and U (r) is the potential 
one [10]. 

The if- function should be linear in order to be constant 
when the coordinates and velocities are being changed. 
It is always possible to represent such function as E = 
T + U by means of scale transformation and usage of 
the necessary system of coordinates. Thus, the sum of 
kinetic and potential energies of the system in a non- 
homogeneous space should be constant. 

Let us take elementary particle with mass m and ve- 
locity v. The kinetic energy corresponding to the particle 
will be T(v 2 ) = mv 2 /2, and the potential energy -U{r), 
so E = mv 2 /2 + U(r) = const. In this case from equality 
E = follows that: 

v(mv + dU/dr) = (1) 

The eq. (1) is a balance equation of the kinetic and 
potential energies. The first term of the equation de- 
termines the change of energy caused by inertness of a 
particle, and the second term is the change of energy 
caused by heterogeneity of space. It is obvious, that the 
equation (1) is carried out, if the condition takes place: 

mv = -dU/dr (2) 

It is Newton equation (NE). This equation is determining 
the connection of acceleration of the particle with the 
external force. The right hand side of eq. (2) is the active 



force. The left hand side is the inertial force [9, 10]. The 
particle moves along the gradient of a potential function. 
The work of forces on the closed line in a potential field is 
equal to zero. The dynamics of the particle is reversible. 

Let us take a system, which consists of N potentially 
interacting elements; the mass of each element is equal to 
1. The force acting on each element is equal to the sum of 
the forces from another elements and the external force. 
The force between elements is central and depends on the 
distance between them. 

The energy of the system is equal to the sum of ki- 

N 

nctic energies of elements - Tjy — £ mv 2 /2, poten- 

i=l 

tial energies in the field of external forces - U~N env , 
and the potential energy of their interaction t/jv(rjj) = 

N-l N 

2~2 2~2 Uij(rij), where = n — rj - is a distance 

i=l j=i+l 

between elements i and j. So, E = En + JJ env = 
T N + U N + U env = const. 

The time derivative of the energy will be as follows: 

N 

= (3) 

i=l 
N 

Where Fi = mhi + £ Fij + Ff nv is effective force for i 
particle; U mv = £ Vi F?™- Fy = dU N /dr l3 ; F?™{n) = 

i=l 

dU env /dn. 

Theeq. (3) can be rewritten as: E — J2 v iFi = 0- This 

i=i 

equality can be treated as orthogonality of the vector of 
effective forces with respect to the vector of velocities 
of elements of the system. If there are no restrictions 
imposed on the V{ directions, the requirement Fi = is 
satisfied [9]. Then from eq. (3) we obtain: 

N 

mv, = -Y^ViFi - Fr v (4) 
»=i 

It is NE for the system's elements in non-homogeneous 
space. From it follows, that the motion of an clement 
of system is determined by the force equal to the sum 
of vectors of forces, acting from all other particles and 
external force [11, 12]. 

Let us obtain an equation of motion of the system as a 
whole in an external field. In variables in which energy of 
system is represented in the form of the sum of energy of 
its motion and the internal energy, the external force act- 
ing on system also is represented us a two corresponding 
forces. For it show, we shall take advantage of equal- 

N N-l N 

ity: T N = £ - m/(2N){V 2 + £ E «?•} 

i— 1 i— 1 

N 

(a), where Vn — Rn — ^/N £ fi -are velocities of the 

»=i 
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center of mass (CM); Rn - are coordinates of the CM; 
v ij = r%j ■ We will write the energy of the system in such 
a way: E N = Tft + E% s , E% s = Tff s + U N . Then the 
eq. (3) can be written as follows: 

N 

Tjv + E% s = - v i F i nv (5) 

i=l 

Where T*f = M N V N V N ;M N = mN; E% s = f^ s + 
UT= E Viiimvij/N + Fij). 

i=l j=i+l 

Let us represent velocity of the motion of elements 
of the system as the sum of velocities of their motion 
with respect to the CM of the system - fj, and veloc- 
ity of the CM of itself - Vn , i.e. Vi — Vi + Vn- Us- 
ing these variables we will have the following: Tn = 



N N N 

£ mv t 2 /2 = m/(2N)V N + mV N E ^ + E mi%/2. As 

i— 1 i— 1 i— 1 

N N 

£ Vi = 0, then T N = m/{2N)V N + £ mvf/2. There- 

i=l i=l 
AT AT-1 AT 

fore £ mvf/2 = 1/(2N) £ £ u?.. Thus the kinetic 

i=l i=l j'=i+l 



energy of the relative motion of particles of the system 
equals the sum of kinetic energies of the particles' motion 
with respect to the CM. 

Let us take into account that = fjj = fj — fj , where 
fi.fj - are coordinates of the elements with respect to 
the system's CM. In this case we can write: [/jv(rjj) = 

N-l N 

UN{fi ) = U N (fi). That is E E % i? y( r y) = 

i=l j=i+l 

N N 

J2viFi(fi), where F { = dU N /dfi = ^ dU N /dr iy 

i=l j^ij=l 

By means of generalization of corresponding equalities 
for the kinetic and potential energies of the system, we 
will obtain from eq. (5): 

AT 

V N M N V N + Y mvifa + F(f)i) = 

i=l 

N 

= -V N F env -J2viF? nv (R,ri) (6) 

i=l 

AT 

Here J F ient ' = ^ Ff~ nv (R,r t ), R - is a coordinate of 

i=l 

CM. 

The eq. (6) determines the balance of energy of the 
system in non-homogeneous space. The first term in the 
right hand side corresponds to change of kinetic energy of 
motion of system as the whole. The second term deter- 
mines the change of internal energy. If the external forces 
are absent the eq. (6) breaks up on independent equa- 
tions: one of them is the equation of motion of the CM, 
and the others are the equations of motion of particles. 

Let us take into account that F env = F env (R + n), 
and suppose that R 3> fj. Then it is possible to ex- 
pand the force F env in a series using a small parameter, 



fi/R. Keeping the terms up to first-order of infinitesimal, 
we will have: Ff nv = Ff nv \ R + {VFf nv )\ R h = Ff™ + 

N N 

(\7Fl^ v )ri. Taking into account that E ^» = E ^* = 

i=l i=l 

N 

and £ F$? v = NFf™ = F ent \ we can set from (6): 

i=l 

N 

V N {M N V N ) + mvi{ii + F(f)j) w 

i=l 

AT 
i=l 

In the eq. (7) the force Fq"" is potential and depends 
on R. It is determines the change of kinetic energy of 
system as the whole. The second term in the right hand 
side is depending on coordinates of particles and their ve- 
locities in relative to the CM. It is determines the change 
of internal energy of system. The force corresponding to 
this term is non-potential and do not change the system's 
momentum as the whole. 

Thus, dynamics of system in an external field is de- 
termined by transformation of two types of energy: the 
energy of motion of system as the whole and its internal 
energy. For each of this type of energy the force de- 
termining its change is corresponds. Energy of motion 
of system is determined by macroparameters - velocities 
and coordinates of the CM. The change of internal energy 
is determined by macroparameters and microparameters. 
Thus, the system is similar to the structured particle. Its 
dynamics is determined by the forces changing its inter- 
nal energy and energy of its motion as the whole. 

Let us note, that the eq. (5) can be obtained, directly 
basing on the NE for elements. For this purpose we shall 
multiply the eq. (4) on the corresponding velocity. After 
summation the obtained equations for all particles we 
shall have the eq. (5) (if we have summarized the eq. (4) 
without multiplying it on velocity in this case the internal 
forces in the second term of the eq. (5) will be lost [13]). 
It is confirms validity of the equation (5). 

III. THE EQUATION OF INTERACTION OF 
TWO SUBSYSTEMS 

The noncquilibrium system can be submitted by a set 
of IES [8] . The dynamics of a set of IES is determined by 
the eq. (6) if the external forces replace on the forces of 
interaction of IES. Thus, to find the equation of motion 
of IES, it is necessary to determine collective forces of 
their interaction. 

Let us the system consists of two IES - L and K. We 
take all elements to be identical and have the same weight 
1, and L to be a number of elements in L - ILS, K - 
is a number of elements in K -ILS, i.e. L + K = N, 

L K 

Vl = 1/i E v i an d Vk = 1/K v i ~ are ILS's velocities 

i=l i=l 
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with respect to the CM of system. The velocity of the 
system's CM we take equal to zero, i.e. LVl + KVk = 0. 

We can represent the energy of the system as En — 
El + E K + U mt = const, where E L and E K are 
the ILS, and U mt - is the energy of their interac- 
tion. According to the eq. (6), the energy of each 
ILS can be represented as E L = Tf + E\ ns , E K 



K 



rptr + E ins^ 

T*T + Ep s , where T\ r = M L Vg/2, T% = M K V*/2, 



Ml = mL, Mk — mK. E ms - is the internal energy 
of a ILS. The E ms consists of the kinetic energy of mo- 
tion of the elements with respect to the CM of IES - 
T ms and their potential energy - U ms , i.e. E ms = 

L-l L 

E E Ui L j L (ri L j L ), 

*i=l 3L=iL+l 

K-l K 



T ins + jjms^ where u%na 



j Tins 
U K 



= E E U iKjK (r iKjK ). The energy U 



is 



IK =1 JK=iK+l 

K L 

determined as U mt = E E Uj L j K {rj L j K ). Indexes 

]K = l ]L = 1 

jk,jL, iK,iL determine belonging of the elements to cor- 
responding ILS. In equilibrium we have: T tr = 0. Hence, 
if the system aspirates to equilibrium, then T tr energy 
for each ILS will be transformed into the internal energy 
of IES. 

We have obtained the equations of dynamics of L and 
K of IES in the following way. Let us differentiate energy 
of system on time. In order to find the equation for L - 
IES, at the left hand side of obtained equality we have 
kept only those terms which determine the change of the 
kinetic and potential energies of interaction of elements of 
L - IES. We replaced all other terms in the right hand side 
and combined the groups of terms in such a way when 
each group contained of the terms with identical veloc- 
ities. In accordance with NE (see eq. (5), the groups 
which contain terms with velocities of the elements from 
K- IES are equal to zero. As a result the right hand 
side of the equation will contain only the terms which 
determine the interaction of the elements L-IES with the 
elements if -IES. The equation for if-IES can be obtains 
in the same way. Then we execute replacement of vari- 
ables: Vi — Vi + V and take into account equality (a). As 
a result we will have [7, 14]: 



L-l L 

v l m l v l +j2 E 

»L = 1 jx,=»L + l 

+F iLJL }} 



ILJL 



+ 



L 



(8) 



K-l K 



V K M K V K +J2 ( v 



IK3K I 



K 



+ 



*jr = l 3K =*jf + 1 

+F iKjK ]} = ^ K + V K ^ (9) 
Here R K = (l/K) £ n K ; R L = (1/L) £ r it ; * = 

i K = i »L = 1 

E E^\ $ L = E v iL Fi[; $ K - £ v iK Ff- K . The 

»L = 1 *L=1 *Jf = l 



terms: F%(R K ,r iL ) = £ F itjjt and Ff K {R L ,r iK ) = 



E ^ 



are forces between the corresponding particle 



of one IES and all particles of the other IES. The work 
of these forces determines the change of energy of IES. 

The eqs. (8, 9) are the EQS. They are describe energy 
exchange between IES. Independent variables EQS are 
macroparameters - coordinates and velocities of motion 
IES, and also microparameters - coordinates and veloc- 
ities of elements. So, EQS connects among themselves 
two types of the description: at a macrolevel and at a 
microlevel. The description at a macrolevel determines 
dynamics IES as the whole, and at a microlevel deter- 
mines dynamics of elements IES. 

The force, determines motion of IES as the whole. 
This force is the sum of the potential forces acting on 
elements of one IES at the side of another IES. 

The forces which determined by the terms <&l and 
transformed of the motion energy of IES to inter- 
nal energy as a result of chaotic motion of elements of 
one IES in a field of the forces of another IES. It is 
non-potential force which cannot be expressed as a gra- 
dient from any scalar function. These forces is equiv- 
alents to dissipative forces. It can be shown with the 
help of some transformations of the eqs. (8, 9). For 

L 

this purpose we take into account that £ mVi [ Vi L = 

L-l L . L 

(V^) E E mv *LjLVt L j L , Ul = E f ^l^ l = 

iL=ljL=iL+l »L=1 

L-l L L 

E E EiLjLVi L j L , F lL = dU L /df iL ■ Then 

EQS can be rewritten so: 



M L V L = -tf - ol l V l 



M k V k = - a K V K 



(10) 



(11) 



where a L = -(E L ns + <5> L )/Vl a K = -{E™ s - 
*x)/V£, E™ = E v iL (mi iL 

»L=1 

K 



+ Fir.), Ef s 



v. 



IK 



{mv iK +F lK ). 



c = l 



Here " a^" , " a/f" are coefficients determining efficiency 
of transformation of energy of relative motion the IES 
into internal energy. They are equivalent to the friction 
coefficients. Thus the role of friction is reduced to redis- 
tribution of energy of motion IES between their elements. 

If the noncquilibrium system is submitted by a set of 
IES, the state of the system can be determined by the 
point in the phase space which consists from 6R — 1 of in- 
dependent variables, where R is a number of IES. In this 
space the role of a elementary particle is carried out by 
the IES. We will call this space as S-space to distinguish 
it from usual phase space with 67V — 1 of independent 
variables. We see that the S'-space, unlike usual phase 
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space, is compressed, though total energy of all elements 
is preserved. The rate of compression of S-space is deter- 
mined by velocity of transformation of energy of relative 
motion the IES into their internal energy The volume 
of compression of S'-space is determined by energy of the 
IES motion. 

Thus, EQS determines the IES dynamics as a result of 
transformation of IES interaction energy into two types of 
energy: internal energy of IES and energy of its motion. 
The forces which are carrying out such transformation 
break up on potential and non-potential parts. Potential 
forces determine the change of the velocity of IES. The 
non-potential force determines the change of their inter- 
nal energy. The work of non-potential forces connected 
with the chaotic motion of elements for one IES in a field 
of forces of another IES. 

IV. THE LAGRANGE, HAMILTON AND 
LIOUVILLE EQUATIONS FOR IES 

Let us explain how Lagrange, Hamilton and Liouville 
equations for the noncquilibrium system which submitted 
by a set of IES can be obtained [14]. The canonical types 
of these equations for elements are following from the 
integral principle of Hamilton [9]. In turn the integral 
principle of Hamilton follows from differential principle of 
D'Alambcrt. D'Alambcrt equation is constructed on the 
basis of NE. In accordance with D'Alambert principle: 
"the virtual work of the effective forces which includes the 
inertial and active forces is equal to zero for all reversible 
virtual displacements of elements compatible with the 
given restrictions" [9]. But if the system consist of a 
set of IES this principle will be as follows: the sum of 
works of all forces of interaction the IES at their virtual 
displacements compatible to restrictions on dynamics is 
equal to zero. 

If the change of the internal energy of the IES can be 
neglected, the work on their motion will be determined 
only by the potential part of interaction forces of IES. In 
this case instead of EQS we can use NE and accept the 
IES as elementary particle. As a result we come to the 
well known canonical equations of a classical mechanics. 

If neglecting by the change of internal energy is impos- 
sible, DAlambert equation should be written down on 
the basis of EQS which take into account the transfor- 
mation of the IES motion energy into internal energies in 
the result of the work of the non-potential part of forces. 
As a result we will obtain Lagrange, Hamilton and Li- 
ouville equations for IES. Let us explain briefly a way 
of these equations of obtaining. Firstly, the DAlambert 
equation on the basis of EQS can be obtained. Basing on 
it we obtain Lagrange equation. After that the Hamilton 
and Liouville equations are deduced. 

Let us take a system consisting of N elements which 
can be represented by a set of IES. Required Lagrange 
equation for noncquilibrium system looks like [14]: 

d/dt{d^/dv n ) - d%/dr n = -F n (12) 



Here 3- is a Lagrange function for the system, n — 
1, 2, 3...N -is a number of the particle. 

At absence of internal degrees of freedom in IES, the 
work of non-potential forces is equal to zero. In this case 
the right hand side of eq. (3.1) is equal to zero and the 
equation becomes canonical [10-12]. 

The Lagrange equations for two L and K of IES can 
be written as: 

d/dt(d<s L /dv iL ) - d$ L /dr iL = -F« (13) 



d/dt(d$ K /dv iK ) - dQ K /dr iK = -Ft K (14) 

Here 3^ and ?$k - are Lagrange functions for the IES. 

The Hamilton's equations for non-equilibrium systems 
can be written as: 

dH/dr n = -p n - F n (15) 



dH/d Pn = v n (16) 

The non-potential part of force of interaction IES to- 
gether with potential force determines the right hand side 
of the eq. (15). 

The Liouville equation for non-equilibrium systems can 
be written as [6, 7]: 

df/dt = fdF n /d Pn (17) 

Her /-is a distribution function for the system's particles. 

The right hand side of eq. (17) is not equal to zero 
as forces between IES depend on velocities. It means 
that S'-space for the noncquilibrium case is compressed. 
When relative motions of IES will disappear, the right 
hand side of eq. (17) will be equal to zero. Thus the 
descriptions in S-space and in the usual phase space for 
equilibrium systems are similar. I.e., at absence of IES 
motion energy, the dissipative processes do not exist. 

V. EQS AND THERMODYNAMICS 

Here we would like to show, how basing on the EQS 
it is possible to come from the classical mechanics to 
thermodynamics. 

Let us take the motionless noncquilibrium system con- 
sisting from "i?" of IES. Each of IES consists from enough 
plenty of elements. Let us, dE is a work which done 
above the system. In thermodynamics it is term as inter- 
nal energy of a system (do not confuse E with the E ms 
- internal energy of IES). The dE is determined by the 
basic equation of thermodynamics as: dE — dQ — PdY 
[8]. Here, according to common terminology, E is energy 
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of a system; Q is thermal energy; P is pressure; Y is 
volume. 

As well as the basic equation of thermodynamics, EQS 
also is differential of two types of energy. According to 
the EQS the volume dE is redistributed inside of the 
system so, that one part of its goes on change of energy 
of relative motion of 1ES, and another part changes their 
internal energy. 

The first term in the left hand side of EQS is a change 
of kinetic energy of motion of a IES as the whole, dT tr . 
This term corresponds to the term PdY. Really [8], 
dT tr = VdV = VVdt = Vdr = PdY. 

If the potential energy is a homogeneous function of a 
second order of the radius- vectors, then as it follows from 
the virial theorem [10], we have: E ins = 2T ins = 2U ms . 
The line denotes the time average. Let us consider 
the system near to equilibrium. The average energy 
of each element is E ins = E lns /N = nT^ ns where N is 
a number of elements. As the increasing of the inter- 
nal energy is determined by the volume dQ, then we 
will have:' dQ w T ins [dE ins /T^ ns ] ~ T^ ns [dv /v ], where 
v is the average velocity of an element, and dvo is its 
change. For the system in the closed volume we have: 
dv /v ~ dT/T, where T is the phase volume of a sys- 
tem, dT will increase due to increasing of the system's 
energy on the value, dQ. By keeping the terms of the 
first order we get: dQ « T* ns dT/T = T™ s d\nT. By def- 
inition rflnT = dS ms , where S ms is a entropy [8]. So, 
near equilibrium we have dE lns = dQ « TQ ns dS lns . 

If we break IES into subsystems, these subsystems will 
not have the relative motion. Therefore the entropy in- 
creasing, AS, for noncquilibrium system is completely 
determined by the energy T tr passing into E ms . There- 
fore AS can be determined by the formula [7, 14]: 

1 = 1 k=l J s 

Here E L is the kinetic energy of L-IES; Nl is the num- 
ber of elements in L-IES; L = 1, 2, 3...R; R is the number 
of IES: s is a number of the external element which in- 
teraction with element k belonging to the L-IES; F^ s is 
a force, acted on fc-element; Vk -is a velocity of the k- 
element. 

Thus in agreement with eq. (18), the entropy is de- 
termined by the energy of the relative motion of IES, 
transformed in an internal energy as a result of work of 
non-potential part of forces between IES. 

To obtain equation for the entropy production wc 
take into account that: AS = AQ/T. Thus dS/dt = 
[dE ms /dt]/(kE ms , where k is a coefficient. It is possi- 
ble to express this formula through the work of forces 
of interaction IES. Let us Eq is a full system's energy, 
E^q is a beginning energy of relative motion of L-IES. 
In according with eqs. (10,11) the rate of increasing of 
the internal system's energy is equal to: £ — 
The internal energy of a system is equal to: E ms — 



E - J2l=i E l, where T,l=i e l is a sum of IES en- 
ergy of relative motion. But X)l=i e< l = Co — J C(t)dt, 
where Co = ^2l—i e To- Then entropy production for 
the system, g pro d = dS/dt, can be write: g pro d — 
D/(l -D + J* D(t)dt), where D = (/E , D = (o/E. 

The energy of T tr characterizes the rate of non- 
cquilibrium system while E ms characterizes its degree of 
equilibrium. To be a stationary state of the system the 
loss of energy E tr is necessary to compensate by the in- 
flow of external energy. It is possible as a result of contact 
of system with thermostat or due to outflow of radiation. 
Then the stationary state of noncquilibrium system is 
characterized by the formula: g pro d = \q~ \ — \ Q+ \ , where 
p_ is entropy outflow, g + -is inflow of entropy, g pro d is 
entropy production, determined by the formula (18). 

Thus, the EQS leads to the consent between classical 
mechanics and thermodynamics. It is possible to explain 
by that the EQS equation determines dynamics of the 
structured particles taking into account the change of 
their internal energy 



VI. CONCLUSION 

Creation of methods of studying of dynamics of the 
open nonequilibrium systems can be carried to a task: 
how to determine the properties and laws of a systems 
dynamics, knowing properties and laws of dynamics of 
its elements, at the external restrictions set on system. 
Within the frame of existing formalisms of a classical 
mechanics its solution collided with the contradiction be- 
tween classical mechanics and thermodynamics. Accord- 
ing to obtained here results the nature of this contradic- 
tion consists in the following. 

Evolution is impossible for systems of elementary par- 
ticles as processes determining it are caused by the dis- 
sipative forces. But such opportunity appears for the 
structural particles possessing an internal energy. The 
change of internal energy causing evolution is related to 
work of non-potential part of the interaction force of par- 
ticles. In many cases: near to equilibrium, in the linear 
approach of the theory of perturbation, etc., the neglect- 
ing of this part of force is admissible. But, it is obvi- 
ous, that at studying evolutionary processes its account 
is necessary. 

Evolution of noncquilibrium system is determined by 
the intrasystem energy streams. Our key idea is that 
these energy streams can be found with the help of EQS 
if the system is submitted by a set of IES. Dynamics of 
IES is determined by the internal energy. Only owing to 
occurrence of this parameter within the frame of classical 
mechanics the determination of entropy is appearing as 
the rate describing increase of internal energy of the IES. 

Independent variables in EQS are macroparameters - 
coordinates and velocities of IES motion and also mi- 
croparameters - coordinates and velocities of elements. 
I.e., EQS connects two types of the description: at 
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a macrolevel and at a microlevel. This description is 
achieved by representation of IES energy by the sum of 
two types of energy: energy of their motion and inter- 
nal energy. The change of each of them due to the IES 
interaction are determined by the corresponding type of 
collective force between of IES. 

The change of the motion energy of IES is determined 
by the work of potential forces. The change of the in- 
ternal energy of IES is determined by the work of non- 
potential forces as a result of chaotic motion of particles. 
The process of transforming of internal energy into the 
energy of IES motion is forbidden by the law of preser- 
vation of momentum. Really, the IES momentum cannot 
be changed by internal motion of its elements. Therefore 
the work on the closed contour for IES is not equal to 
zero due to increasing of the internal energy of IES. 

The Lagrange, Hamilton and Liouville equations on 
the basis of EQS are deduced. These equations describe 
dynamics of the noncquilibrium system in the compress- 
ible iS-space where points is corresponds to the momen- 
tum and coordinates of the IES. 

The following mechanism of irreversibility can be pro- 
posed: the energy of relative motion of IES is transformed 
into their internal energy as a result of the work of the 
non-potential part of interaction force between of IES. 
The system equilibrates, when relative motion of IES dis- 
appears. 

The explanation of the First law of thermodynamics is 
based on the fact that the work of subsystems' interaction 
forces changes both the energy of their motion and their 
internal energy. The explanation of the Second law of 
thermodynamics is based on the condition of irreversible 
transformation of the subsystems' relative motion energy 
into their internal energy. 



The dynamics of nonequilibrium systems is determined 
by other form of symmetry, than dynamics of the equi- 
librium systems. Really, both for elementary particles, 
and for equilibrium systems the dynamics is reversible 
because internal energy are absent. But for nonequilib- 
rium systems these types of symmetry are broken. It is 
due to the work of the non-potential forces of IES interac- 
tion that change the internal energy. Thus it is possible 
to assume, that the broken of symmetry (for example, 
observable broken CP-symmetry in disintegrations (3 - 
meson [15]), can testify that products of disintegration 
are the structured particles with time of a life commen- 
surable with the period of supervision. 

If we take into account that all particles are struc- 
tured we come to the conclusion about existence of infi- 
nite hierarchical sequence of the microstructures enclosed 
in macrostructures. In turn, a macrostructure are a mi- 
crostructure for the following hierarchical step. The oc- 
currence and existence of such hierarchical sequence in 
the nature is provided with hierarchy of fundamental 
forces of interaction. As a rule, the force of interaction of 
the bottom hierarchical level there are more than forces 
of the top level. Indeed, nuclear forces more than elec- 
tromagnetic forces and electromagnetic forces more than 
the gravity force. We come to the infinite divisibility of 
a matter because of impossibility of occurrence of ele- 
mentary particles. This conclusion has the experimental 
confirmation [15, 16]. 

Thus the offered approach allows us to expand a scope 
of formalism of classical mechanics on the open nonequi- 
librium systems and eliminate contradictions between 
classical mechanics and thermodynamics. 
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